DPNU-94-32 
KEK Preprint 94-101 
TIT-HPE-94-07 
TUAT-HEP 94-3 
INS-REP-1062 
NWU-HEP 94-03 
OCU-HEP 94-5 
KOBE-HEP 94-04 
PU-94-684 



Measurement of the forward- backward asymmetries 
for charm- and bottom-quark pair productions at 
(y^)=58GeV with electron tagging Q 



TOPAZ Collaboration 
E.Nakano"' |] R■Enomoto^ M.Iwasaki^ KAbe", TAbe'^, IAdachi^ KAdachi^ MAoki'^, 
MAoki", SAwa^ R.Belusevic^ K.Emi^ H.Fujii^ K.Fujii^T.Fujii-^, J.Fujimoto^ K.Fujita^ 
N.Fujiwara^ H.Hayashii^ B.Howell'', N.Iida^ H.Ikeda^ R.Itoll^ Y.Inoue^, H.Iwasaki^ 
K-Kaneyuki"^, R.Kajikawa", S.Kato*, S.Kawabata^ H.Kichimi^ M.Kobayashi^ D.Koltick'^, 
I.Levine'*, S.Minami'^, K.Miyabayashi", A.Miyamoto'', K.Muramatsu'^, K.Nagai-', K.Nakabayashi", 
O.Nitoh% S.Noguchi^ A.Oclii^, F.Oclliai^ N.Ohislli^ Y.01lnislli^ Y.Olisliima^ H.Okuno*, 
T.Okusawa^,T.Shinoliara^, A.Sugiyama", S.Suzuki", S.Suzuki'', K.Takahaslii^, T.Takaliaslii^, 
T.Tanimori'^, T.Tauclii'', Y.Teramoto^, N.Toomi'^, T.Tsukamoto'', O.Tsumura*^, S.Uno'', 
T.Watanabe'', Y.Watanabe'', A.Yamaguchi'^, A.Yamamoto'',and M.Yamauchi'' 



*to be published in Physics Letters B 

tinternet address;NAKANOE@KEKVAX.KEK.JP 



1 



(a) Department of Physics, Nagoya University, Nagoya 464, Japan 

(b ) KEK, National Laboratory for High Energy Physics, Ibaraki-ken 305, Japan 

(c) Department of Physics, Nara Women's University, Nara 630, Japan 

(d) Department of Physics, Tokyo Institute of Technology, Tokyo 152, Japan 

(e) Department of Applied Physics, Tokyo Univ. of Agriculture and Technology, Tokyo 184, Japan 

(f) Department of Physics, University of Tokyo, Tokyo 113, Japan 

(g) Department of Physics, Osaka City University, Osaka 558, Japan 

(h) Department of Physics, Purdue University, West Lafayette, IN 47907, USA 

(i) Institute for Nuclear Study, University of Tokyo, Tanashi, Tokyo 188, Japan 

(j) The Graduate School of Science and Technology, Kobe University, Kobe 657, Japan 
(k) Faculty of Liberal Arts, Tezukayama University, Nara 631, Japan 



Abstract 

We have measured, with electron tagging, the forward-backward asymmetries 
of charm- and bottom-quark pair productions at (^/i)=58.01GeV, based on 23,783 
hadronic events selected from a data sample of 197pb^^ taken with the TOPAZ 
detector at TRISTAN. The measured forward-backward asymmetries are A^^ = 
-0.49±0.20(stat.)±0.08(sys.) and A^.^ = -0M±0M{stat.)±0A3{sys.), which 
are consistent with the standard model predictions. 
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1 Introduction 



The differential cross section for fermion pair productions, e^e —>■ ff, can be 
written in tlie following form in the massless limit: 



dajf- 3 



-aff(l + cos^ 9 H — Apn cos( 



dcosO 8 3 

where a^j and Ap^ are the total cross section and the forward-backward charge asym- 
metry, respectively, while 6 is the polar-angle of the final-state fermion / with respect 
to the direction of the initial-state electron. 
In the standard model|]l|, Ap^ is given by 
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with 

s 

^ ~ 16 sin^ Ow cos2 Owis - M|o + iMzoT%o) ' 
where Ve{vf) and ae{af) are the vector and axial-vector coupling constants of the elec- 
tron (final-state fermion) to the boson, Q/ is the charge of the final-state fermion, 
and Mzo and F^o are the mass and the total width of the boson, respectively. The 
standard model predicts 

^;^ = 2(//-2g/sin2 Ow) 

and 

aj = 2li, 

where // is the third component of the weak isospin of the fermion (/) and 9]^ is the 
Weinberg angle. 

The above formula tells us that Ap^ attains to its maximum in the TRISTAN energy 
region and that its measurement there is sensitive to and therefore to the structure of 



the multiplet to which the fermion belongs. The measurement of the forward-backward 
charge asymmetry thus provides a good test of the standard modeL The predicted 
asymmetries for charm- and bottom-quark pair productions are 

A'pj, = -0.47, 
= -0.59 

at yi=58.01GeV for M^o = 91.1888GeV/c2, r^o = 2.4974GeV, and sin^ = 0.2321 
0. The A^p^ is, however, reduced by the B-B mixing, whose probabihty x is given by 
X = RdXd + RsXs^ where and Xi are the fraction and the mixing probabihty of Bd 
or Bs, respectively. Their measured values are 

Xd = 0.16, = 

and 

0.3 

Xs = 0.53, i?, = — , 

so that the expected Ap^ value becomes —0.43. 

In this energy region, only the data from the TRISTAN experiments are available 
and the TRISTAN average values of A^^g |, |] and A^^^, g |, |, |ig have been 



A'p^ = -0.56 ±0.09, 
= -0.59 ±0.09. 

The Ap^ is consistent with the standard model predictions and the A^p^ is slightly 
deviated from it. 

The TRISTAN data include our previous measurement |^ of the forward-backward 
charge asymmetry for the charm-quark pair production through both exclusive and 
inclusive reconstructions of D*^ n^D'^: Ap^ = — 0.49lQj3(stat.) ± 0.06{sys.), con- 
sistent with the standard model prediction. 



In order to improve the statistical accuracy, we have carried out another measure- 
ments using electron tagging to be described in this paper. It should be noted that 
this measurement is completely independent of the above D* analysis and that, in our 
energy region, there is only one previous A'^p^ measurement reported, which are used 
lepton tagging at this y/s p. 

2 The TOPAZ detector 

The main components of the TOPAZ detector[Tl[] include a time projection cham- 
ber (TPC) and a barrel lead-glass calorimeter (BCL) which were essential to electron 
identification. 

Combining the dE/dx information from the TPC and the E/P information from the 
BCL, we were able to select electrons in hadronic final states with high purity and high 
efficiency over a broad momentum range. 

Instead of getting into details of these detectors, we summarize their performance 
here. The momentum resolution of the TPC has been measured to be 

-ft 



through e"'"e~ fi^fi~ and cosmic fi^ events|ll], while its dE/dx resolution was deter- 
mined to be 4.6% by a study of minimum ionizing pions. The energy resolution of the 
BCL can, on the other hand, be expressed as 



E \\^E{GeV) 

3 Analysis 
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3.1 Electron selection 

This analysis is based on 23,783 hadronic events. The selection method was 
described in Ref[lT2|. This data sample corresponded to an integrated luminosity of 
197pb"^ and was taken at an averaged center-of-mass energy (y^)=58.01GeV. 

In search of electron track candidates, we first selected good charged tracks from 
the hadronic events, using the following selection criteria defining a good track. 

1. The closest approach to the interaction point (i?) had to be less than 1.0cm in 
the X-Y plane (perpendicular to the beam axis) and that in the Z direction (Z) 
to be less than 4.0cm. 

2. the absolute value of the cosine of the polar angle had to be between 0.02 and 
0.83, 

3. the transverse momentum (P^) with respect to the beam axis had to be greater 
than 0.15GeV/c, 

4. the number of degrees of freedom (N.D.F.) in the track fitting had to be greater 
than 3, and 

5. there had to be more than 30 hit wires for dE/dx calculation (65%-truncated 
mean) out of 114 wires maximum. 

Each good track was extrapolated to the BCL to look for its corresponding BCL 
cluster and to test E/P. The clustering of the energy deposits in the BCL was carried 
out iteratively by merging a counter to its neighboring counter if its energy was smaller 
than that of the neighboring counter. For each of so formed clusters, we calculated its 
energy as the sum over counters and its position as the energy-weighted mean of the 
counter positions. 
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Figure 1 : (a) The distribution of the minimum distance between the extrapolated track position on 
the BCL and the cluster position, (b) The distribution of the cluster width of a shower with respect 
to the track-extrapolated position. The solid and dashed histograms are for the electron-enhanced 
sample defined in the text and for all tracks, respectively. 

Comparing the extrapolated track position and the BCL cluster positions, we looked 
for the cluster that is the closest to the track. The closest distances between tracks 
and clusters are histogrammed in Fig.^(a) for all tracks (dashed) and for an electron- 
enhanced sample (solid) selected by requiring 0.75<E/P<1.25 and 5.5<dE/dx<7.5(keV/cm) 
We accepted those tracks which had a distance less than 5.0 cm. 

We also calculated the energy-weighted r.m.s. of counter positions in the cluster 
(cluster width) with respect to the matched track position, whose distribution is shown 
in Fig.|l|-(b). The selected tracks were further required to have a cluster width between 
1.0 and 10.0cm. 

Since background tracks were predominantly in the low-momentum region, we im- 
posed an additional momentum cut P > 0.8GeV/c. 

Finally, the dE/dx information from the TPC was used to complete our electron 
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selection: the for electron hypothesis {xD had to be xl < 3.0 (N.D.F. = 1). The 
tracks rejected by this cut (dE/dx- rejected hadrons) are to be used in the background 
estimation. 



3.2 Rejection of e^e pair background 

We rejected electrons apparently coming from 7 conversions or Dalitz decays as 
follows. Secondary vertices were reconstructed for all unlike-sign pairs of the TPC 
tracks. When the two tracks in a pair did not intersect in the X-Y plane, we required 
the pair to have a distance at the closest point less than 4.0cm in the X-Y plane, and 
2.0cm in the Z direction. The pair also had to have a deflection angle in the X-Y plane 
less than 5.0 degrees between its momentum vector and the flight direction from the 
interaction point. If the invariant mass of the pair was less than 50MeV/c^, the tracks 
were rejected. 

When the tracks intersected in the X-Y plane, we selected from the two intersections 
the one which gave the smaller deflection angle as the secondary vertex. We then 
required the pair to have a distance less than 2.0cm in the Z direction and to have a 
deflection angle less than 5.0 degrees. If the invariant mass was less than 200MeV/c^, 
the tracks were rejected. 

In this way, most of the electrons coming from 7 conversions or Dalitz decays were 
removed. Nevertheless, there still remained a significant number of electrons from 7 
conversions or Dalitz decays, which were estimated through Monte-Carlo simulation |T^: 
the Monte-Carlo simulation gave us the ratio of the number of all the reconstructed 
conversion tracks to that of the remaining tracks. 

Using this ratio, we estimated the number of remaining pair-conversion tracks from 
the actual number of all the reconstructed pairs in the experimental data on a bin by bin 
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basis and subtracted them from the electron candidates. By doing this, we can reduce 
the systematic errors due to the error of the material thickness in the detector simula- 
tion program. The remaining electrons from 7 conversions or Dalitz decays estimated 
through the Monte-Carlo simulation are 311.7 ± 17.7 events. The E/P distribution for 
these electron candidates is shown in Fig.|[ 

3.3 Hadron background 

The hadron background was estimated using the dE/dx- rejected hadrons in the 
experimental data and is shown in Fig.0 as the dashed histogram, whose normalization 
factor was calculated so as to equalize the entries in the side-band (E/P = 0.0 - 0.64) 
for the electron candidates and the background sample. The normalized background 
was subtracted from the electron candidates and the remaining electrons in the region 
0.72<E/P<2.00 were counted. This method had been checked out through the Monte- 
Carlo simulation. 

The estimated number of electrons from primary charm, b-to-c cascade, and direct 
bottom decays are 151.5 ± 12.3, 66.5 ± 8.2, and 131.9 ± 11.5 events, where the selection 
efficiencies of 10.4%, 11.3%, and 24.1%, respectively. 

3.4 Charm- and bottom-quark sample 

We calculated the transverse momentum (Pt) of each electron candidate track with 
respect to the axis of jets reconstructed using the invariant-mass algorithm 0]. For the 
jet reconstruction, we used charged tracks with momenta greater than 0.2GeV/c and all 
of the neutral clusters, which were clusters with a distance to the closest track greater 
than 5.0 cm. 

We studied the angular resolution of the jet axis with respect to the primary quark 
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Figure 2: The E/P distribution of the selected electron samples. The solid histogram shows the 
experimental data. The dashed histogram is for the estimated background by dE/dx-rejected hadrons. 

direction through the Monte-Carlo simulation. Since the resolution of the thrust axis 
was about 9 degrees, while that of the reconstructed jets as above was about 6 degrees, 
we used the reconstructed jets as the primary quark directions. 

The Pt of an electron from a charm-quark is expected to be lower than that of a 
bottom-quark in general. To enhance charm or bottom contents, therefore, we divided 
the sample at Pt = 0.8GeV/c into two classes: low-P^ (charm-enhanced) and high-P^ 
(bottom-enhanced) samples. 

3.5 Monte Carlo simulation 



We generated Monte-Carlo events using the JETSET6.3 generator [|T3[ with M^o 



91.173GeV/c2,r^o = 2.487GeV, and ?>m^ Ow = 0.2325[|. For light quark events, its 
parameters were tuned by a multi-parameter fit of hadronic event shapes[|l^]. For 
heavy quark events, we adjusted the parameters for fragmentation function, so as to 
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match other experiments^, to be a=0.8 and b=0.2. The B-B mixing effect is included 
in our Monte-Carlo simulation. Using this Monte-Carlo simulation, we estimated the 
acceptance and radiative correction factors to be used later in the following subsection. 



3.6 Fitting procedure 

Figs.|^(a) and -(b) show the momentum and the Pt distributions of electrons, 
respectively. The —Q cos 6 distributions of the high-P^ and the low-Py electrons are 
shown in Figs.^(a) and -(b), respectively, where Q is the charge of the electron and 6 is 
the polar angle of the jet axis with respect to the electron beam axis. The points with 
error bars in the figures are experimental data, while the histograms are the best-fit 
results obtained by fitting Fig.|^(a), Figs.^(a), and -(b) simultaneously. The used fit 
function to the —Qcos6 distribution is 



N^F— ■ 2Br(c ^e)F( 



<?<? 



+ N;;^P—-2Br{b 



+ N'^P^-2Br{c 



e)F{A^ 



FBJi 



■ cf 



e)(l + x)F( 



X 



1 + x 



A 



FB)i 



■Cf 



with X = Br{b ccs) set to be 16%||T^, where, Ng^^ is the number of hadronic events, 
and F[A-'p'^)i are integrals of the following formula over i-th cos 6' bin; 

3 



'l + cos^e + -A 



FB 



cost 



Cf^, cf'', and Cf"" are the Monte-Carlo-determined correction factors for the i-th bin of 
electrons from prompt charm, prompt bottom, and b-to-c cascade decays, respectively, 
which were described in the previous subsection and they are listed in Table |I[ 
The fit determined the branching fractions of c e and 6 — » e to be 



Brie 



0.131 ±0.015 
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Table 1: List of correction factors. 




Figure 3: (a) The momentum distribution of e='='s. (b) The Pt distribution of e='='s. The points 
are background-subtracted data. The open, hatched, and double-hatched histograms are the best-fit 
results for the contributions from c e+X, b ^ e+X, and b e~X, respectively. 
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Figure 4: The -Qcos0 distributions of e='='s (a)in the low- Pt region {Pt < 0.8GeV/c) and (b)in the 
high-Pr region (Pt > 0.8GeV/c). The points are experimental data. The open, hatched, double- 
hatched, triple-hatched, and filled histograms are the best-fit results for the contributions from c 
e^X, b c ^ e^X, c — > e~X, 7 conversion, and Dalitz decay, respectively. 

Br{h-^e) = 0.109 ±0.025 

and and A^p^ to be 

A'pj^ = -0.49 ±0.20 

= -0.64 ±0.35. 

3.7 Systematic errors 

We checked various systematic-error sources. The estimated systematic errors are 
summarized in Table ^ and the varied parameter values used for the estimation are listed 
in Table 0. The dependence on the selection of good tracks was checked by changing 
the cut values on R,Z,Pt, and momentum (P). The dE/dx dependence was tested by 
changing the xl cut. The error caused by the shower shape parameters was estimated by 
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error source 


ABr{c~*e) 


^FB 


ABr{6->e) 
Br{b^e) 


^^'fb 
^fb 


track selection 


7.4% 


3.3% 


6.4% 


17.7% 


dE/dx 


1.0% 


7.5% 


1.3% 


3.2% 


cluster selection 


5.7% 


6.8% 


10.5% 


0.6% 


background 


2.2% 


8.2% 


1.7% 


6.7% 


7 and Dalitz rejection 


1.5% 


1.0% 


1.6% 


1.5% 


Pt cut 


3.0% 


8.3% 


8.1% 


7.2% 


total 


10.2% 


15.8% 


14.9% 


20.6% 



Table 2: Summary of systematic errors. 



changing the distance cut and the cluster width cut. We checked the effect of the back- 
ground estimation by changing the region for background normalization and electron 
counting: the first set used 0.0<E/P<0.64 and 0.72<E/P<2.00 for normalization and 
counting, respectively, while the second one used 0.40<E/P<0.60 and 1.6<E/P<2.00 
for normalization and 0.72<E/P<1.28 for counting. The systematic error due to the 
estimation of the 7 conversions and Dalitz decays was checked by changing the cuts 
for the 7 conversion and Dalitz decay rejection. The QCD correction was estimated by 
comparing the JETSET 6.3 P.S. and JETSET 6.3 qq options to be -0.02 and -0.01 for 
and A^p^, respectively. Since they were small, we neglected them. The errors from 
the Pt cut is also given in Table ^. The errors due to uncertainty of B-B mixing were 
also checked and they were negligibly small(~ 0.2%). The overall systematic errors 
were obtained by adding them in quadrature. 

3.8 Results and discussion 

Our results are 

Br{c e) = 0.131 ± 0.015(stal) ± 0.013(s?/s.) 
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error source 


parameter nominal value 


varied value 


track selection 


R 


1 .0cm 


1.2cm 




Z 


4.0cm 


3.0cm 




Pt 


0.15GeV/c 


0.20GeV/c 




P 


0.8GeV/c 


l.OGeV/c 




xl 


3.0 


9 7 


cluster selection 


distance 


5.0cm 


7.5cm 




width 


1.0 - 10.0cm 


2.0 - 12.5cm 


7 and Dalitz rejection 


X-Y distance 


4.0cm 


5.0cm 




Z distance 


2.0cm 


2.5cm 




deflection angle 


5 degrees 


7 degrees 


Pt cut 




0.80GeV/c 


0.75GeV/c 


B-B mixing 


X 


0.139 


0.0 


Table 3: 


Summary of varied parameter values. 




Br{b ^ e) = 


0.109 ±0.025(stat.) 


±0.016(s?/s.) 




^FB ■ 


= -0.49±0.20(stat. 


) ±0.08(s?/s.) 






= -0.64±0.35(stot.) 


1 ±0.13(sys.). 





The branching fractions are consistent with the previous measurements, and the forward- 
backward asymmetries are consistent with the standard model predictions as well as 
our previous measurements |, |^. The obtained forward-backward asymmetry for 
charm-quark pair production is plotted in Fig. ^ together with other experimental data 

ill, 1,0- 

The combined result with our previous D* analysis is 

A'pj^ = -0.49 ±0.12, 

from which we obtain the charm-quark's axial-vector coupling constant(a^) to the 
boson to be = 1.07lo!3i5 being consistent with the standard model. 
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4 Conclusion 

We have measured the forward-backward asymmetries of charm- and bottom- 
quark pair productions via e~^e~ annihilations through an inclusive electron analysis 
at (-yi)=58.01GeV. The number of hadronic events used for this analysis is 23,783, 
corresponding to an integrated luminosity of 197pb^^. 

The measured forward-backward asymmetries are Ap^ = —0.49 ± 0.20{stat.) ± 
0.08{sys.) and A^p^ = —0.64 ± 0.35(stat.) ± 0.13(sys.), consistent with the standard 
model prediction and our previous measurements. Combining the results from our 
previous D*"^ measurement, we obtained Ap^ = —0.49 ± 0.12 and = l.OTloisf. 
The obtained branching fractions are Br{c ^ e) = 13.1 ± 1.5{stat.) ± 1.3{sys.)% and 
Br{b ^ e) = 10.9 ± 2.5{stat.) ± 1.6{sys.)%, which are in good agreement with the 
previously measured values of 9.6 ± 0.9%||T^ and 10.8 ± 0.5%[1^, |T^, respectively. 



16 



5 Acknowledgement 

The authors thank Dr. M. Sakuda for discussions on the analysis. They also thank 
the TRISTAN accelerator staff for the successful operation of TRISTAN. The authors 
appreciate all of the engineers and technicians at KEK as well as the collaborating 
institutions: H. Inoue, N. Kimura, K. Shiino, M. Tanaka, K. Tsukada, N. Ujiie, and H. 
Yamaoka. 



17 



References 

[1] S. Weinberg, Phys. Rev. Lett. 19 (1967) 1264; 

A. Salam, in: Elementary particle theory; relativistic group and analyticity, Nobel 
Symp. No. 8, ed. N. Svartholm (Almqvist and Wiksells, Stockholm, 1968) p. 361; 
S. L. Glashow, Nucl. Phys. 22 (1961) 1579. 

[2] D. Schaile, in: 27th Intl. Conf. on High Energy Physics (Glasgow, 1994) 

[3] Particle Data Group, Phys. Rev. D45, No.ll, Part II (1992). 

[4] E. Nakano et al., Phys. Lett. B 314 (1993) 471. 

[5] K. Abe et al., Phys. Lett. B 313 (1993) 288. 

[6] F. Hinode et al., Phys. Lett. B 313 (1993) 245. 

[7] F. Liu et al., KEK Preprint 93-147 submitted to Phys. Rev. D 

[8] A. Shimonaka et al., Phys. Lett. B 268 (1991) 457. 

[9] K. Nagai et al., Phys. Lett. B 278 (1992) 506. 

[10] M. Shirakata et al., Phys. Lett. B 278(1992) 499. 

[11] R. Enomoto et al., Nucl. lustrum. Methods A 269 (1988) 507; 
A. Imanishi et al, Nucl. lustrum. Methods A 269 (1988) 513; 
T. Kamae et al., Nucl. lustrum. Methods A 252 (1986) 423; 
A. Yamamoto et al, Jpn. J. Appl. Phys. Lett. 25 (1986) L440; 
S. Kawabata et al., Nucl. lustrum. Methods A 270 (1988) 11; 
J. Fujimoto et al, Nucl. Instrum. Methods A 256 (1987) 449; 



18 



S. Noguchi et al., Nucl. Instrum. Methods A 271 (1988) 404; 
K.Fujii et al, Nucl. Instrum. Methods A 264 (1988) 297. 

[12] 1. Adachi et al., Phys. Lett. B 255 (1991) 613. 

[13] T. Sjostrand, Comput. Phys. Commun. 39 (1986) 347; 

T. Sjostrand and M. Bengtsson, Comput. Phys. Commun. 43 (1987) 367. 

[14] J. Leveille, in: Proc. second Moriond Workshop (Moriond, 1982), eds. J. Tran 
Thanh Van and L. Montanet, P.191 

[15] M. Althoff et al., Phys. Lett. B 126 (1983) 493; 
W. Braunschweig et al., Z. Phys. C 44 (1989) 365; 
W. Bartel et al., Phys. Lett. B 146 (1984) 121; 
F. Ould-Saada et al, Z. Phys. C 44 (1989) 567; 
H. Aihara ct al., Phys. Rev. D 34 (1986) 1945; 
M. Derrick et al., Phys. Lett. B 146 (1984) 261; 
P. Baringer et al, Phys. Lett. B 206 (1988) 551; 
H. Aihara et al, Z. Phys. C 27 (1985) 39; 
H. Aihara et al., Phys. Rev. D 31 (1985) 2719; 
M. Althoff et al., Phys. Lett. B 146 (1984) 443; 
H. J. Behrend et al., Z. Phys. C 47 (1990) 333; 
E. Elsen et al., Z. Phys. C 46 (1990) 349; 
B. Adeva et al., Phys. Rep. 109 (1984) 131. 

[16] H. Albrecht et al., Phys. Lett. B 249 (1990) 359. 

[17] S. Henderson et al., Phys. Rev. D 45 (1992) 2212. 



19 



